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Abstract 

Background: Arabinogalactan proteins (AGPs) are cell wall proteoglycans that have been shown to be important 
for pollen development. An Arabidopsis double null mutant for two pollen-specific AGPs {ogpG agpll) showed 
reduced pollen tube growth and compromised response to germination cues in vivo. A microarray experiment was 
performed on ogp6 ogpll pollen tubes to search for genetic interactions in the context of pollen tube growth. 
A yeast two-hybrid experiment for AGP6 and AGPl 1 was also designed. 

Results: The lack of two specific AGPs induced a meaningful shift in the gene expression profile. In fact, a high 
number of genes showed altered expression levels, strengthening the case that AGP6 and AGPl 1 are involved in 
complex phenomena. The expression levels of calcium- and signaling-related genes were found to be altered, 
supporting the known roles of the respective proteins in pollen tube growth. Although the precise nature of the 
proposed interactions needs further investigation, the putative involvement of AGPs in signaling cascades through 
calmodulin and protein degradation via ubiquitin was indicated. The expression of stress-, as well as signaling- related, 
genes was also changed; a correlation that may result from the recognized similarities between signaling pathways in 
both defense and pollen tube growth. 

The results of yeast two-hybrid experiments lent further support to these signaling pathways and revealed putative 
AGP6 and AGPl 1 interactors implicated in recycling of cell membrane components via endocytosis, through 
clathrin-mediated endosomes and multivesicular bodies. 

Conclusions: The data presented suggest the involvement of AGP6 and AGPl 1 in multiple signaling pathways, in 
particular those involved in developmental processes such as endocytosis-mediated plasma membrane remodeling 
during Arabidopsis pollen development. This highlights the importance of endosomal trafficking pathways which are 
rapidly emerging as fundamental regulators of the wall physiology. 
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Background 

Pollen-pistil interaction is initiated when the male gameto- 
phyte is transferred from the anther to the stigma. The 
pollen grain then starts to hydrate and germinate, forming 
a pollen tube that grows through the carpels internal 
tissues to deliver its two sperm cells into the embryo sac. 
Pollen tubes elongate through the extracellular matrix of 
the pistil tissues, extending by an actin-myosin-based tip- 
growth mechanism that transports vesicles loaded with 
new cell wall material to the extending apex and finally, 
when arriving at the embryo sac, enter into the synergid 
cell. Fundamentally during this type of localised, tip- 
focused growth, the cytoplasm is polarized directing 
secretory events to the tip. Internal gradients and trans- 
membrane ion fluxes, notably of calcium ions, are another 
key feature of pollen tube growth [1]. 

Arabidopsis pollen transcriptome analysis showed the 
expression of a unique subset of genes relative to the 
sporophytic tissues [2-6]. Furthermore, changes in gene 
expression patterns occur during the development of the 
male gametophyte from a microspore to a mature tricellu- 
lar pollen grain [7]. More recently, genome-wide expres- 
sion profiling of pollen tubes grown in vitro identified 
another set of genes that are expressed in the pollen tube 
but not in pollen, which suggests de novo mRNA synthesis 
in the growing pollen tube [5]. The gene expression pro- 
files of in vitro- and semi in v/vo-grown Arabidopsis pollen 
tubes have also been characterized and were found to dif- 
fer; this lead to the discovery of a specific sub-set of genes 
that are activated by potentiation of the pollen tube by the 
pistil [6]. 

Among the genes, or gene families that may be involved 
in the transition from a sporophytic to a gametophytic- 
developmental program are those which code for arabino- 
galactan proteins (AGP). AGPs constitute a large family of 
cell wall proteoglycans that are found on the plasma mem- 
brane, in the cell wall, in the apoplastic space, and in secre- 
tions. They have also been found in detergent-resistant 
membranes isolated from Arabidopsis suggesting their 
presence in lipid rafts [8]. Key distinguishing features of 
AGPs are: (1) the carbohydrate, usually branched type II 
arabino-3,6-galactan, 0-linked to Hyp residues of the pro- 
tein backbone, which constitutes 90-98% (w/w) of the 
molecular mass, (2) the protein backbone typically rich in 
the dipeptide motifs Ala-Hyp, Ser-Hyp, Thr-Hyp, Val-Pro, 
Gly-Pro and Thr-Pro, (3) a glycosylphosphatidylinositol 
(GPI) membrane anchor, predicted to be present on most, 
but not all, AGPs and (4) the ability to bind to a class of 
synthetic chemical dyes, known as Yariv reagents [9], 
which are not only useful for detection, quantification and 
isolation of AGPs, but also for fiinctional studies [10]. 

The identification of Arabidopsis pollen-specific AGP 
genes has been used to scrutinize phenotypic changes in 
the respective null mutants. We previously identified 



two male gametophyte AGP genes {AGP6 and AGP 11) 
which showed functional redundancy [11]. An agp6 
agpll double null mutant was subsequently obtained. 
The double null mutant exhibited segregation distortion, 
as assessed by the number of aborted pollen grains, sug- 
gesting that the gametophyte generation was affected. 
The strong reduction in pollen germination and pollen 
tube growth rate, together with premature ectopic ger- 
mination of pollen (whilst it was still in the anther) [12], 
prompted us to further analyze this mutant, in order to 
gain insight into the mode of action of AGPs. We 
decided to examine the transcriptome of the agp6 agpll 
double null mutant pollen tubes, using the Affymetrix 
ATHl Genome Array. 

Here we report the identification of 1022 genes whose 
expression in the double null mutant pollen tubes was 
shown to be either reduced or elevated, when compared 
to wild type pollen tubes. These genes can be used as 
starting points to dissect the gene regulatory networks 
in which AGPs are involved during pollen tube growth. 
In parallel, we performed yeast two-hybrid experiments 
to identify interactors of AGP6 and AGPll, and to pro- 
vide evidence for the biological functions of these AGPs. 

Results 

An Arabidopsis line simultaneously null for two pollen- 
specific AGP genes, AGP6 and AGPll, has been charac- 
terized in our laboratory and found to display notable 
phenotypic alterations, namely partially aborted pollen 
grains, reduced germination potential and precocious ger- 
mination inside the anthers [11,12]. Furthermore, under 
the conditions used in this study pollen tubes of the 
double mutant were statistically shorter (83 [im) than 
those of the wild type (100 (im). 

To further characterize the agp6 agpll mutant line a 
differential microarray screen was carried out to identify 
genes with modified expression in the 8 h-grown pollen 
tubes of agp6 agpll compared to wild type pollen tubes. 



General data analysis 

The number of expressed genes in the wild type pollen 
tubes (as indicated by a minimum of 2 out of 3 Present 
calls in the triplicate GeneChip experiments) was found to 
be 6886, which is approximately 28% of the total number 
of genes represented in the GeneChip Arabidopsis ATHl 
Genome Array. This value is in close agreement with 
those published for microarray studies of Arabidopsis 
pollen tubes using the same type of microarray chips and 
experimental set up, i.e., in vitro-gvown pollen tubes har- 
vested 4-6 h post-germination [5,6]. Microarray data was 
published in GEO (http://www.ncbi.nlm.nih.gov/geo/info/ 
linking.html) with the record number GSE40861. 
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Functional classification of wild type pollen tube 
transcriptomes 

The transcriptome of wild type pollen tubes (6886 genes) 
was categorized into 34 functional groups according to 
the MapMan visualization software (Figure 1). Of the 
genes indicated as present, 28.7% had unknown functions, 
whereas the rest were primarily involved in protein syn- 
thesis (18.8%), RNA transcription and processing (10.6%), 
signaling (5.0%), transport (5.0%), and cell organization 
and cytoskeleton (4.9%). This functional distribution is 
consistent with the recruitment of cell resources for pollen 
tube growth, and is in agreement with other pollen and 
pollen tube microarray studies [4,5]. 

Arabidopsis agp6 agpll differentially expressed genes 

The number of genes expressed in the agp6 agpll mutant 
pollen tubes was roughly equal to that of the wild type 
control although only 87.5% of expressed genes were com- 
mon to both data sets (Figure 2). There are two gene popu- 
lations of specific interest; 1) those present only in the 
double mutant and 2) those present only in the wild type 
pollen tubes. However, a statistical analysis of the signifi- 
cantly differentially expressed genes necessarily produced 
a rather different population of genes from that shown 
in Figure 2. Genes were considered to be differentially 



expressed from those they were compared with if the 90% 
lower confidence bound of the fold change between experi- 
ment and baseline was above 1.3, resulting in a median 
False Discovery Rate (FDR) of less than 5%. Observing 
such criterion a population of 1022 differentially expressed 
genes was obtained which accounts for 14.7% of the pollen 
tube transcriptome (Figure 2 and Additional file 1). These 
1022 genes included 155 genes present only in agp6 
agpll, 168 genes present only in wild type, and 699 
genes expressed in both. About 60% of these 699 genes 
were up-regulated in agp6 agpll, with the remaining 
40% being down-regulated. Overall the agp6 agpll mu- 
tation caused an up-regulation of over 500 genes in 
Arabidopsis pollen tubes. 

Validation of gene expression 

RNA levels were independently verified for a collection of 
genes, either by real-time PGR (qPGR) or by conventional 
semi-quantitative RT-PGR. Genes were selected on the 
basis of expression differences between the two microarray 
data sets and of absolute signal values (Table 1). The rela- 
tive expression levels for all genes tested were consistent 
with the results of the microarray experiments, thus con- 
firming its reliability. 



Not assigned 
Protein 
RNA 
Signaling 
Transport 
Cell 

Miscelaneous 
Stress 
DNA 

Development 
Cell wall 
Lipid metabolism 
Hormone metabolism 
Redox 

Amino acid metabolism 
Secondary metabolism 
IVIitochondrial electron transport / ATP synthesis 
Nucleotide metabolism 
Minor CHO metabolism 
PS 

TCA / Org. transformation 
Major CHO metabolism 
Glycolysis 

Co-factor and vitamin metabolism 
Metal handling 
Cl-metabolism 
Biodegradation of xenobiotics 
Tetrapyrrole synthesis 
OPP 

Polyamine metabolism 
N-metabolism 
S-assimilation 
Fermentation 
Gluconeogenesis / Glyoxylate cycle 
Micro RNA, natural antisense etc 
No classification 




Figure 1 Functional classification of wild-type Arabidopsis pollen tube transcriptome. Analysis was performed with Classification 
SuperVi ewer Tool using MapMan source [13]. 
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Figure 2 Comparison of thie wild-type and agp6 agpll pollen tube transcriptomes. Number of genes identified as being Present' in 
wild-type pollen tubes (WT/P) and in ogp6 ogpll pollen tubes (DM/P) are shown on the left, and genes "Present" among the differentially 
expressed genes in wild-type pollen tubes {WT_D\f_P) and agp6 agpll pollen tubes (DM_Dif_P) are shown on the right. Among the whole set 
of 1022 differentially expressed genes, 155 were only present in the double mutant, and 168 genes were only present in the wild-type. The 
remaining pool of 699 genes contained both those which were up- and down-regulated. 



Functional classification of differentially expressed genes 

According to the MapMan classification software, the 
functional classification of the entire double mutant 
transcriptome, does not substantially diverge from the 
functional classification of the wild type transcriptome 
(Figure 1). However, a significant number of genes are 
differentially expressed, i.e. they are either up- or down- 
regulated in response to the absence of two AGP gene 
products. The MapMan functional classification of these 
differentially expressed genes is shown in Figure 3. 

The MapMan functional BINs "Protein", "RNA", 
"Transport", "Signaling", "Cell", "DNA, and "Stress" are 



the groups in which the greatest number of genes was 
altered, and together account for almost half of the dif- 
ferentially expressed genes (46%). 

A sub-cluster of the larger "Protein" cluster that con- 
tains genes with F-box motifs was particularly striking for 
its differential signal intensities and number of genes 
affected, with a representation factor of 1.8 x relative to the 
Arabidopsis genome. Likewise, the sub-clusters "signaling, 
calcium", and "stress.abiotic.heat" show enrichment fac- 
tors of 1.5 X and 2.1 x, respectively, in the differentially 
expressed gene set, as compared to the whole genome. 
The sub-clusters "signaling.receptor.kinases" and "stress. 



Table 1 Summary of confirmatory PCR assays 



Gene 


AGI ID 


LBFC 


FC - qPCR'' RT-PCR^ 


Ubiquitin-conjugating enzyme 20 


Atlg50490 


+2.44 


+7.0 


Proline transporter 1 


At2g39890 


+2.02 


+2.0 


Zinc finger (C3HC4-type RING finger) family protein 


At5g60250 


+ 1.96 


+4.0 


Sucrose phosphate synthase 2 F 


AtSgllllO 


+ 1.70 


+ 1.5 


Expansin B5 


At3g60570 


+ 1.50 


+2.0 


CML42 


At4g20780 


+ 1.48 


+4.5 


Plant cadmium resistance 1 1 


Atlg68610 


-1.99 


-2.0 


Cellulase 3 


Atlg71380 


-1.98 


-1.4 


Phosphoinositide 4-kinase gamma 4 


At2g46500 


-1.75 


-2.2 


Heat shock protein 17.6A 


Atlg53540 


-1.34 


-2.2 


CAP (Pathogenesis-related protein, putative) 


At2g 19970 


+23.34 


+ 


C2 domain-containing protein 


At3g57880 


+ 16.93 


+ 


F-box family protein 


Atlg65760 


+ 11.8 


+ 


Beta glucosidase 36 


Atlg51490 


+4.95 


+ 


CAP (allergen V5) 


At2g 19980 


+4.79 


+ 


Hypothetical protein 


At2g22340 


-141.79 




Pentatricopeptide (PPR) 


At5g28380 


-44.55 




Avirulence-responsive family protein 


At4g09950 


-25.77 




DET3 


Atlg 12840 


-12.71 





^ Lower bound of fold change (microarray experiment). 
^ agp6 agpl 1 to wild-type fold change for qPCR experiments. 
Genes up- or down-regulated In agp6 agpl 1, as assessed by semi-quantitative RT-PCR are indicated with a (+) or a (-) sign, respectively. 
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Figure 3 Classification of differentially expressed genes according to MapMan software. The proportion of genes, up- or down- regulated, 
within a given gene classification type is indicated and the actual number of genes this represents is shown. 



biotic" were not enriched in the differentially expressed 
group of genes but were selected because of the high 
LBFC value of some of its members (Table 2). 

Another group of stress-related genes (CAP genes; 
Table 3) was also notable for their signal intensities and, 
in particular the large increase in the respective signal 
intensities of two of its members compared to wild type 
(Table 3). 

To identify possible gene interactions and/or associa- 
tions we used a Web interface (GeneMANIA) which is 
designed to analyze gene lists and prioritize genes for 
functional assays [15]. Extracting the differentially 
expressed genes from the list of 500 genes with the high- 
est signal intensity, a gene network was produced that 
highlighted a "response to heat" gene cluster (Figure 4). 

AGP6 And AGP1 1 interactors detected by yeast-two 
hybrid assays 

To experimentally detect putative interactors for AGP6 
and AGPll we used the core domains of these two pro- 
teins as baits to screen a normalized Arabidopsis expres- 
sion library. A whole plant cDNA library was used as it 
was thought to be the best choice for identifying AGP 
interactors. The library was generated from pools of 
cDNAs obtained by mRNA extracted from several 
organs at different developmental stages including devel- 
oping inflorescences, developing siliques, and mature 
flowers just before, and at, anthesis. The library was nor- 
malised in order to achieve a better representation of 
low abundant messengers. 



For the AGP6 library mating plates. 111 colonies grew 
on selective media. Forty-eight among them were selected 
randomly for sequencing and they revealed 37 unique pos- 
sible candidates. The ones that had an absent call in our 
array experiment as well as in other publicly available 
pollen and pollen tube arrays [5-7] were excluded, result- 
ing in a set of 22 likely candidates for interaction (Table 4). 
For the AGPll library mating plates only 11 colonies grew 
on selective media and all were sequenced. Those resulted 
in 5 likely candidates for interaction (Table 4). 

Of the putative interactors uncovered, four were selected 
to confirm interaction with AGP6 and one with AGPll: 
CaLB domain family protein (Atlg70810), MAPK9 
(At3gl8040), UBQIO (At4g05320) and FBR12 (Atlg26630), 
and AP2C1 (At2G30020), respectively. All co-transformed 
colonies with the AGP and each candidate grew on all se- 
lective media, while each negative control did not, thus 
confirming the interaction (Figure 5). 

Discussion 

Sexual plant reproduction relies on the transfer of pollen 
to the stigma, where it hydrates and germinates into a 
pollen tube that extends through the pistil tissues to its 
target, the embryo sac. The agp6 agpll double null mu- 
tant, besides showing a pollen tube growth defect also 
shows an early "inside the anther" germination pheno- 
type. These deficiencies, provide further evidence for the 
role of AGPs in cell wall deposition and growth, and po- 
tentially represent a defect in receiving, or responding 
to, germination cues. 
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Table 2 agp6 agpll differentially expressed genes among the MapMan BINs [29.5.11.4.3.2: protein, degradation, 
ubiquitin. E3. SCF. F-box], [30.2: signaling.receptor kinases], [30.3: signaling.calcium], [20.1:stress.biotic], 
and [20.2.1: stress.abiotic.heat] 



MapMan BIN AGI locus 


LBFC^ 


Gene annotation'^ 


29.5.1 1.4.3.2: protein, degradation, ubiquitin. E3. SCF. F-box 


Up-regulated 


Atlg65760 


11.8 


Protein of unknown function (DUF295) 


At2g27520 


6.63 


F-box and associated interaction domains-containing protein 


At5g38390 


5.43 


F-box/RNI-like superfamily protein 


At3g52030 


3.92 


F-box family protein with WD40/Y\^N repeat domain 


At2g24080 


3.52 


Protein of unknown function (DUF295) 


At5g27750 


3.19 


F-box/FBD-like domains containing protein 


At3g49520 


2.79 


F-box and associated interaction domains-containing protein 


Atlg70360 


2.65 


F-box family protein 


At3g 17710 


2.40 


F-box and associated interaction domains-containing protein 


At2g 17020 


2.25 


F-box/RNI-like superfamily protein 


At3g 13820 


2.08 


F-box and associated interaction domains-containing protein 


At2g03580 


1.99 


F-box family protein-related 


At5g03000 


1.75 


Galactose oxidase/kelch repeat superfamily protein 


Atlg06630 


1.71 


F-box/RNI-like superfamily protein 


At3g50080 


1.61 


VFB2, VIER F-box proteine 2 


At5g60610 


1.55 


F-box/RNI-like superfamily protein 


At3g 18980 


1.50 


EIN2 targeting proteinl (ETPl) 


At3g59000 


1.44 


F-box/RNI-like superfamily protein 


At2g07120 


1.38 


F-box associated ubiquitination effector family protein 


At5g49610 


1.36 


F-box family protein 


Atlg55590 


1.35 


RNI-like superfamily protein 


Atlg23780 


1.34 


F-box family protein 


Down-regulated 


At5g43190 


-7.32 


Galactose oxidase/kelch repeat superfamily protein 


Atlg68050 


-4.00 


"Flavin-binding, kelch repeat, F-box 1" 


Atlg 12490 


-3.55 


F-box associated ubiquitination effector family protein 


Atlg55000 


-1.99 


Peptidoglycan-binding LysM domain-containing protein 


Atlg 78840 


-1.96 


F-box/RNI-like/FBD-like domains-containing protein 


At2g25490 


-1.58 


EIN3-binding F-box protein 1 


At4g04690 


-1.52 


F-box and associated interaction domains-containing protein 


Atlg70590 


-1.44 


F-box family protein 


30.2: signaling, receptor kinases 


Up-regulated 


At4g23280 


2.27 


Cysteine-rich RLK (RECEPTOR-like protein kinase) 20 


At3g21970 


2.22 


Domain of unknown function (DUF26) 


At5g40380 


1.33 


Cysteine-rich RLK (RECEPTOR-like protein kinase) 42 


At4g20790 


1.76 


Leucine-rich repeat protein kinase family protein 


At3g24550 


1.61 


Proline extensin-like receptor kinase 1 


At4g28670 


1.41 


Protein kinase family protein with domain of unknown function (DUF26) 


At5g23170 


1.35 


Protein kinase superfamily protein 


Down-regulated 


At4g20530 


-7.70 


Receptor-like protein kinase-related 


At5g59650 


-3.61 


Leucine-rich repeat protein kinase family protein 
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Table 2 agp6 agpll differentially expressed genes among the MapMan BINs [29.5.11.4.3.2: protein, degradation, 
ubiquitin. E3. SCF. F-box], [30.2: signaling.receptor kinases], [30.3: signaling.calcium], [20.1 :stress.biotic], 
and [20.2.1 : stress.abiotic.heat] (Continued) 





At4g02420 


-2.87 


Concanavalin A-like lectin protein kinase family protein 




Atlgll280 


-2.01 


S-locus lectin protein kinase family protein 




Atlg61380 


-1.73 


S-domain-1 29 


30.3: signaling, calcium 


Up-regulated 




Atlg 18890 


2.66 


Calcium-dependent protein kinase 1 




At5g 19360 


2.30 


Calcium-dependent protein kinase 34 




At5g62390 


1.57 


BCL-2-associated athanogene 7 




At4g20780 


1.48 


Calmodulin-like 42 Calcium; sensor involved in trichome branching 




Atlg05150 


1.44 


Calcium-binding tetratricopeptide family protein 




At5g42380 


1.40 


Calmodulin-like 37 




Atlg51960 


1.38 


IQ-domain 27 




At5g47100 


1.31 


Calcineurin B-like protein 9 


Down-regulated 




Atlg 18530 


-1.89 


EF hand calcium-binding protein family 




At5g55990 


-1.88 


Calcineurin B-like protein 2 




At2g33990 


-1.68 


IQ-domain 9 




Atlg62480 


-1.54 


Vacuolar calcium-binding protein-related 




At2g41410 


-1.48 


Calcium-binding EF-hand family protein 




At4g34150 


-1.33 


Calcium-dependent lipid-binding (CaLB domain) family protein 


20.1: stress, biotic 


Up-regulated 




At2g 19970 


23.34 


CAP (Cysteine-rich secretory proteins, Antigen 5, and Pathogenesis-related 1 protein) superfamily protein 




At3g02840 


2.56 


ARM repeat superfamily protein 




At4g 16920 


2.52 


Disease resistance protein (TIR-NBS-LRR class) family 




Atlg32210 


1.78 


Defender against death (DAD family) protein 




At2g35520 


1.77 


Defender against death (DAD family) protein 


Down-regulated 




At4g09950 


-25.77 


P-loop containing nucleoside triphosphate hydrolases superfamily protein 




At3g44400 


-6.01 


Disease resistance protein (TIR-NBS-LRR class) family 




At4g 16860 


-4.48 


Disease resistance protein (TIR-NBS-LRR class) family 




Atlg75040 


-2.63 


Pathogenesis-related gene 5 




At3g 18690 


-1.88 


MAP kinase substrate 1 




Atlg 14530 


-1.82 


TOM THREE HOMOLOG 1 (THHl); Protein of unknown function DUF1084 




Atlg58170 


-1.59 


Disease resistance-responsive (dirigent-like protein) family protein 




At2g25240 


-1.41 


Serine protease inhibitor (SERPIN) family protein 


20.2.1: stress, abiotic, heat 


Up-regulated 




Atlg28210 


2.40 


DNAJ heat shock family protein 




At3g22530 


1.97 


Unknown protein 




Atlg72070 


1.50 


Chaperone DnaJ-domain superfamily protein 


Down-regulated 




At2g29500 


-4.93 


HSP20-like chaperones superfamily protein 




At3g46230 


-1.92 


Heat shock protein 17.4 




At5g 12030 


-1.84 


Heat shock protein 17.6A 




Atlg79920 


-1.84 


Heat shock protein 70 (Hsp 70) family protein 
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Table 2 agp6 agpll differentially expressed genes among the MapMan BINs [29.5.11.4.3.2: protein, degradation, 
ubiquitin. E3. SCF. F-box], [30.2: signaling.receptor kinases], [30.3: signaling.calcium], [20.1 :stress.biotic], 
and [20.2.1 : stress.abiotic.heat] (Continued) 



Atlg59860 


-1.78 


HSP20-like chaperones superfamily protein 


At5g22060 


-1.65 


DNAJ homologue 2 


Atlg71000 


-1.49 


Chaperone DnaJ-domain superfamily protein 


At2g03020 


-1.49 


Heat sliocl< protein HSP20/alplia crystallin family 


At4g 13830 


-1.43 


DNAJ-like 20 


At5g02500 


-1.41 


Heat shock cognate protein 70-1 


Atlg53540 


-1.34 


HSP20-like chaperones superfamily protein 



^ agp6 agpll to wild type lower bound of fold change. 
^ TAIR annotations [14]. 



AGPs are a class of molecules whose study has been 
particularly challenging. More than 90% of the molecular 
mass is carbohydrate (predominantly consisting of arabin- 
ose and galactose residues), and the role of the polypeptide 
chain, if any, is unknown. Moreover, how the sugar con- 
tent and quality varies between different gene products, 
and whether different gene products are identically glyco- 
sylated is equally poorly characterised. A direct conse- 
quence of this is that techniques and approaches used to 
study proteins often do not apply to AGPs. The massive 
sugar component which surrounds the protein core is 
likely to have a biological role but its biological interac- 
tions and mechanism of action have so far eluded charac- 
terisation. Phenotypic differences between wild type and 
double mutant pollen tubes can only be attributed to the 
lack of both AGP6 and AGPll. It is impossible to deter- 
mine whether the differential gene expression is due to a 
direct effect of the lack of AGPs on the expression of these 
genes, or whether the altered expression profile is a sec- 
ondary and unspecific consequence of the reduction in 
growth rate. However, an expression-based study should 
improve our understanding of AGPs and their functionality 
by identifying genes that are altered in response to the lack 
of specific AGPs. This type of investigation is likely to be 
particularly informative when applied to a process known 
to involve AGPs, in a plant structure, the pollen tube, 
whose transcriptome is less than a third of the genome. 

This study yielded a number of candidate genes that 
may directly or indirectly interact with AGPs and have a 
role in pollen tube growth. However, a comparative ana- 
lysis of the present data with the expression profile of 
Arabidopsis culture cells subjected to Yariv phenylglyco- 
side reagent, which specifically binds to and precipitates 
AGPs [44-48], was considered to be a valuable aid to the 
identification of genes which may be directly affected by 
the lack of AGP6 and AGPll. 

Microarray screen 

The overall difference obtained for expressed genes was 
substantial, with 1022 differentially expressed genes in 



the double null mutant, compared with the wild type 
control. These 1022 differentially expressed genes could 
be further classified into three subgroups: 155 genes that 
were expressed "de novo" in the double mutant, 168 
genes that were absent in the mutant but present in the 
control, and a group of 699 genes whose expression 
levels were found to be altered but that did not affect 
the present/absent status. 

The gene clustering using MapMan software shows 
that the cluster "Protein" contains the largest fraction of 
differentially expressed genes. "Protein degradation" is 
undoubtedly the most represented group, in which the 
category of ubiquitin-mediated protein degradation is 
evident, with several of those genes belonging to the 
F-Box protein family. The F-box is a protein motif of ap- 
proximately 50 amino acids that function as a site for 
protein-protein interaction [49]. Arabidopsis alone was 
reported to contain approximately 700 F-Box genes [50]. 
By comparison, only 20, 27, and 69 F-Box proteins are 
encoded by Saccharomyces cerevisiae, Drosophila mela- 
nogaster, and human genomes, respectively [51]. Such an 
apparent massive investment has also been found to be 
consistent with genetic studies linking specific plant F- 
Box proteins to numerous processes, including hormone 
perception and signaling, stress protection, chromatin 
remodeling, homeostasis, circadian rhythms, self-incom- 
patibility, and defense against pathogens [52]. In Arabi- 
dopsis, recent genetic studies provided unequivocal 
evidence for the role of this protein degradation system 
during gametogenesis, which involved F-Box 17 [53] and 
a pair of ubiquitin-specific proteases (UBP3/UBP4) [54]. 

The analysis of the data obtained in the present work 
also shows calmodulins (CaMs), in the cluster of signal- 
ing by calcium, among the differentially expressed genes 
with highest LBFC score, namely At5g24880 and 
At4g20780 (CML42), the latter functioning in trichome 
branching and being strongly expressed in pollen as 
shown by CML42::GUS transgenic plants [55]. Calcium 
signals are known to play important roles in plant 
growth and development, including pollen tube growth. 
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Table 3 Summary of expression data for genes belonging to the CAP (Cysteine-rich secretory proteins. Antigen 5, and 
Pathogenesis-related 1 protein) superfamily protein 



AGI locus 


agp6 agp1 1 signal 


Wild type signal 


LBFC^ 


Wild type call'' 


Pollen tube call 


Pollen specificity 


At2g 19970 


299 


7.5 


+ 23.0 


A 


A 


No 


At2g 19980 


5622 


887 


+ 4.8 


P 


P 


Yes 


AtSg 19690 


4785 


4400 


<1.3 


P 


P 


Yes 


AtlgOlSlO 


5490 


4100 


<1.3 


P 


P 


Yes 


At4g25780 


4578 


4036 


<1.3 


P 


P 


Yes 


At3g09590 


1367 


1395 


<1.3 


P 


P 


No 


At5g02730 


672 


686 


<1.3 


P 


P 


No 



^ agp6 agpll to wild type lower bound of fold change. 
A = absent; P = present. 



and CaMs are the most ubiquitous calcium sensors in 
eukaryotes. Popescu et al [56] developed a protein micro- 
array to comprehensively investigate Arabidopsis CaM/ 
CaM-like interactions and identified new targets for these 
calcium sensors, one of which was an F-Box protein. The 
relationship between AGPs and calcium is unsurprising, 
as these proteins are involved in pollen tube cell wall 
growth, and calcium is known to regulate pollen tube 
growth and fertilization [57]. The data obtained may sug- 
gest an interaction of AGPs with calcium through cal- 
modulin. This was also implied using an alternative 
approach to study AGP function which involved the use 



of (p-D-Glc)-3 Yariv reagent that specifically binds to, and 
presumably inactivates, a wide range of AGPs. Guan and 
Nothnagel [44] analyzed the gene expression profile of 
Arabidopsis culture cells subjected to Yariv phenylglyco- 
side reagent for 1 and 10 h in a microarray study. Despite 
the specificities of the biological materials, several differ- 
entially expressed genes were found to be common be- 
tween both experiments (Additional file 2). The values of 
LBFC used in the present work could not be used for dir- 
ect comparison with the data of Guan and Nothnagel [44] . 
However, all values of fold change of agp6 agpll genes 
shown in Additional file 2, have values of LBFC above the 
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Figure 4 Heat shock protein genes are differentially expressed in the agp6 agpll pollen tubes. The internet tool GeneMANIA was used to 
produce a gene association networl<. Tlie input data was all the differentially expressed genes isolated from the 500 genes with highest signal 
intensities in the wild-type microarray data set. Some of those genes (Atlg53540, HSP20-like chaperones superfamily protein; At3g46230, 17.4 kDa 
class I heat shock protein; AtSg 12030, heat shock protein 17.6A; At2g29500, HSP20-like chaperones superfamily protein) belong to a "response to 
heat" functional cluster (pink circles) [15]. 
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Table 4 Interactors of AGP6 and AGP1 1 identified by Yeast two-Hybrid library screening assays 



AvjI IOCUS 


1 Dcr"^ 
Lbrv. 


Gene description'^ 


Expression 
in mature 
pollen*^ 


Expression 
in pollen 
tube^ 


Subcellular localization*^ 


Other 


AGP6 Atlg62750 


-1.66 


Elongation factor G (SCOl) 






Apoplast [16] 




Atlg70810 


-1.02 


Calcium-dependent lipid-binding 
(CaLB domain) family protein 


[5-7] 


[5,6] 






At5g05010 


-1.00 


Clathrin adaptor complexes 
medium subunit family protein 

ilVIUU \d.\di 


[5-7] 


[5,6] 


Plasmodesmata [17], 
cytosol [18], plasma 
membrane [1 9] 


Predicted GPI-attachment 
[20], predicted secretory 
pathway [21], predicted 
beta-barrel TM domain [22] 


AIogUODjU 


— U.yo 


A 1 r-citrate lyase 1 


rv! 
L/J 


LbJ 






Atl 932640 


-0.95 


MYC-related transcriptional 
activator (MYC2) 




[5,6] 




No predicted binding site 
(cacatg), but two similar 
sites (acacaag) in AGP6 
promoter region [23] 


Atlg67090 


-0.93 


Ribulose bisphosphate 
carboxylase small chain lA 

inoLo 1 Aj 


[5-7] 


[5,6] 


Apoplast [16] 


Predicted secretory 
pathway [21] 


At5g55070 


-0.91 


Dihydrolipoamide 
succinyltransferase 


[5], [6] 


[5,6] 


Mitochondrion [24] 




At2g20580 


-0.86 


26S proteasome regulatory 
subunit Nl (RPNIA) 


[6] 


[6] 


Cytosol [18], 
plasmodesmata [17] 


One to five predicted 
alpha-helix TM domains 
[22], ubiquitin-binding 
subunit of proteasome 
complex [25] 


At3g61190 


-0.85 


BON 1 -associated protein 
1 (BAPl) 


[5], [7] 


[5,6] 


Plasma membrane [26] 




At3g09440 


-0.83 


Heat shock protein 70-3 


[5-7] 


[5,6] 


Cytosol [18], plasma 
membrane [27,28], cell wall 
[29], tonoplast [30], Golgi 
[31], apoplast [16] 


Predicted secretory 
pathway [21] 


At5g21105 


-0.77 


L-ascorbate oxidase 


[6] 


[6] 


Plasmodesmata [17], cell 
wall [32] 


Predicted secretory 
pathway [21], predicted 
alpha-helix TM domains 
[22] 


At5g28540 


-0.52 


Luminal-binding protein 1 (BIPl) 


[5-7] 


[5,6] 


Cell wall [29], cytosol [18], 
plasmodesmata [17], 
plasma membrane [27], ER 
[33], tonoplast [30], 
cnioropiast pzj 


Predicted secretory 
pathway [21], one to two 
predicted alpha-helix TM 
domains [22] 


At5g57560 


-0.43 


Xyloglucan 

endotransglucosylase/hydrolase 

ZZ [ 1 L.n4, A 1 nZZ) 




[6] 


Cell wall [24] 


Predicted secretory 
pathway [21], predicted 
alpha-helix TM domain and 
predicted beta-barrel TM 
domain [22] 


At5g 15090 


0.55 


Mitochondrial outer membrane 
protein porin 2 (VDAC3) 


[6] 


[6] 


Plasma membrane [19,28], 
cell wall [29] 


Predicted beta-barrel TM 
domain [22] 


Atl g 16920 


0.71 


Ras-related small GTP-binding 
protein (RABAlb) 


[7] 








At3g 18040 


0.75 


Mitogen-activated protein kinase 
9 (MAPK9) 


[5-7] 


[5,6] 




Predicted GPI-attachment 
[20], predicted alpha-helix 
TM domain [22] 


At4g05320 


0.84 


Polyubiquitin 10 (UBQIO) 


[5-7] 


[5,6] 






At2g20990 


0.87 


Synaptotagmin A {SYTA) 


[7] 


[6] 


Endosome [34], plasma 
membrane [19,28,35], 
ER [32] 


Predicted secretory 
pathway [21], predicted 
alpha-helix TM domain [22] 


Atlg26630 


0.89 


Eukaryotic translation initiation 
factor 5A-2 (FBR12) 


[5-7] 


[5,6] 




Predicted 

GPI-attachment [20] 


At3g 10740 


0.92 




[6], [7] 


[6] 


Apoplast [16], cell wall [36] 
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Table 4 Interactors of AGP6 and AGP1 1 identified by Yeast two-Hybrid library screening assays (Continued) 







Alpha-L-arabinofuranosidase 
1 (ASDl) 








Predicted secretory 
pathway [21], one to two 
predicted alpha-helix 
TM domains [22] 


At5g 11200 


1.22 


DEAD-box ATP-dependent 
RNA helicase 56 


[6] 


[5,6] 


Cell wall [29], 
plasmodesmata [17] 




Atlg31150 


2.50 


Uncharacterized protein 


[7] 






Predicted secretory 
pathway [21], predicted 
alpha-helix TM domain [22] 


AGPll At3g47810 


-1.43 


MAGI Homolog of yeast 
retromer subunit VPS29 


[5-7] 


[5,6] 


Multivesicular body [37], 
endosome [38] 


Mediator of protein 
targeting to vacuole [39] 


Atl 942990 


-0.97 


bZIP transcription factor 
60 {BZIP60) 


[5-7] 


[5,6] 




Predicted GPI-attachment 
[20], predicted alpha-helix 
TM domain [22] 


Atlg32640 


-0.95 


MYC-related transcriptional 
activator (MYC2) 




[5,6] 




Two predicted binding 
sites in AGPl 1 promotor 
region [23] 


At2g 19070 


-0.63 


Spermidine liydroxycinnamoyi 
transferase (SHT) 








Tapetum-specific [40], 
involved in exine formation 
[40,41] 


At2g30020 


0.74 


Putative protein pliospliatase 
2C 25 


[5-7] 


[5,6] 


Nucleus [42,43], 
cytoplasm [43] 





^ agp6 agp11 to wild type lower bound of fold change. 
^ TAIR gene annotations [14]. 
^ in published microarray experiments, 
empirically determined in published references. 



threshold of 13, which was the criterion used in our ana- 
lysis. Genes altered in Yariv-treated cells also include sev- 
eral heat-shock protein genes, which is a functional group 
also found to be affected in agp6 agpll. 

Comparative pistil transcriptome analyses [58,59] identi- 
fied a number of genes encoding proteins potentially 
involved in stress and defense responses. This is particu- 
larly interesting as there is evidence that the molecules 
involved in pollination and stress/defense responses may 
be evolutionary and functionally related. Many authors 
have proposed that certain self-incompatibility mechan- 
isms may have arisen through the modification of pre- 
existing pathogen defense mechanisms [60] . 

The analysis of the stress cluster of the present data 
shows that the genes associated with biotic and abiotic 
types of stress, namely PR and heat shock protein genes, 
had their expression levels significantly shifted. In this 
cluster, genes such as At2gl9970 and At2gl9980 code 
for proteins of the CAP (cysteine-rich, antigen 5 and 
pathogen-related 1 protein) superfamily. Members of the 
CAP superfamily are widespread in living organisms, 
[61], and have been intensely studied in mammalian 
fertilization mechanisms. They are most often secreted 
proteins and are involved in several processes including 
cell adhesion during fertilization. The plant pathogenesis 
proteins of the PR-1 family, which are synthesized dur- 
ing pathogen infection or other stress-related responses, 
belong to the CAP superfamily but the precise functions 
of these proteins are still unresolved. In a recent study 



[62] it was suggested that the glioma pathogenesis - 
related 1 (GLIPRILI), along with other members of the 
CAP superfamily and several other proteins, are involved 
in the binding of sperm to the oocyte complex. Also, in 
a work of microarrays of rice embryo sac cells [63], it 
was revealed that allergen V5/Tpx-l-related proteins, 
which are members of the CAP superfamily, are abun- 
dant in the rice synergid cell. Collectively these findings 
strengthen the possible role of CAP domain-containing 
proteins in cellular adhesion and fertilization. 

The high level of expression of 7 members of the CAP 
superfamily in Arabidopsis pollen tube is relevant given 
their predicted signaling function. Moreover, knowing 
that these proteins are most likely secreted and the fact 
that one of its members (At2gl9970) is only present in 
the agp6 agpll double mutant, may indicate some type 
of relationship between AGPs and CAP processing. 

Application of a different bioinformatic tool, GeneMania, 
a co-expression network was obtained which highlighted 
a stress functional gene cluster. Surprisingly, some of 
the genes correlated with those whose expression was 
found to be altered by treatment with Yariv reagent for 
10 h (At2g29500, At3g46230, Atlg53540 and At5gl2030) 
[44]. Guan and Nothnagel [44] compared the genes that 
were induced by Yariv treatments with genes whose ex- 
pression had been previously been shown to be induced by 
other conditions and concluded that the gene expression 
profile induced by Yariv treatment was similar to that of 
wound response. Knowing that pollen tube growth and 
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Figure 5 Yeast two-hybrid experiments with candidate genes selected from the library screening. YSD media -tryptophane, -leucine, 
-liistidine 1 mM 3-Amino-l,2,4-triazole. Empty AD and BD vectors were used as negative controls. (A) Assays with AGP6 protein core as bait and 
CaLB domain family protein (Atlg70810), MAPK9 (At3gl8040), UBQIO (At4g05320) and FBR12 (Atlg26630) as preys. (B) Assays with AGP]] protein 
core as bait and AP2C1 (At2G30020) as prey. 



discharge and the defense against fungal attack are alike in 
many respects, and that both pollen tubes and fungi exploit 
similar receptor proteins [64,65] we propose a role for 
AGPs in this process where, apparently, pollen tube and 
fungal hyphae activate plant cell responses. 

Yeast two-hybrid experiments 

To identify possible AGP6 and/or AGPll interactors, 
yeast two-hybrid experiments were performed. These 
experiments revealed some interesting partners involved 
in the proteosome-independent roles of ubiquitination 
in signaling and endocytosis. 

Besides the well-known role in proteosome degradation, 
ubiquitin conjugation is also involved in down-regulation 
of membrane receptors, transporters and channels. Ubi- 
quitination of plasma membrane proteins leads to their 
endocytosis into the multivesicular endosome and most 
members of the ubiquitin ligase family responsible for traf- 
ficking of diverse proteins carry an N-terminal calcium- 
dependent lipid/protein C2 domain that specifically binds 
phosphoinositides in yeast [66,67]. 



Endosomes are primarily intracellular sorting organelles, 
and receive proteins and lipids from both the biosynthetic 
and the endocytic pathways. Plasma membrane proteins 
that are internalized by endocytosis are either recycled 
back to the plasma membrane or sorted for degradation, 
which is achieved by intermediate/late endosomes, also 
called multivesicular bodies (MVBs) [68,69]. 

It is known that pollen tubes grow by rapid tip loca- 
lized exocytosis, most probably coordinated with an also 
tip localized endocytosis. Using Arabidopsis and Nicoti- 
ana pollen tubes as models, Zhao et al. [70] showed that 
phosphatidylinositol-4-phosphate 5-kinase 6 (PIP5K6) 
regulates clathrin-dependent endocytosis in pollen tubes. 

The internalization and secretion of arabinogalactan- 
rich glycoproteins through MVBs has been described 
but only microscopically [71,72]. The data obtained in 
the present study suggests a physiological interaction of 
AGPs with members of the endosomal system, and that 
AGP6 and AGPll, two specific and highly expressed pro- 
teins in Arabidopsis pollen and pollen tubes may be 
involved in this signaling pathway that sends and recycles 
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proteins to the extracellular matrix. In the yeast two- 
hybrid screening assays, clathrin adaptor complexes were 
present such as MUG13.13 and MAPK9, SYTA (a trans- 
membrane protein involved in membrane trafficking), 
R ABA lb (a small GTP binding protein also involved in 
intracellular protein trafficking) and MAGI (a homolog of 
the yeast retromer subunit VPS29). These results give fur- 
ther support for the clathrin machinery of receptor intern- 
alization in pollen tubes. 

Although AGPs are presumably acting extracellularly, 
they are targeted to the outer leaflet of the plasma mem- 
brane through vesicles, which originate in the ER and 
Golgi network. During the trafficking into the extracellular 
matrix it is probable that AGPs interact with cytosolic 
proteins such as bZIP60 or FBR/eIF5A, two regulatory 
proteins that were identified in the yeast two-hybrid 
experiments. 

With endocytosis established as an essential plant cell 
function [69] and AGP6 and AGPll supposedly involved 
in this process, we looked further into the microarray 
data. Three Rab GTPases known as specific molecular 
markers of plant endosomes, Ara6 (At3g54840), Ara7 
(At4gl9640) and Rhal (At5g45710, also a member of 
the heat stress transcription factor family) are all down- 
regulated in agp6 agplL RabA2 (At3g46830) and Rab A3 
GTPases define a trans-Golgi endosomal membrane that 
overlays with VHA-al (At2g28520) one of the isoforms 
of a membrane integral V-ATPase specifically localized 
to the TGN [73]. Another subunit of V-ATPases, DET3, 
is one of the most down- regulated genes in the present 
study and known to be involved in plant growth and de- 
velopment [74]. 

A C2 domain is a protein structural domain involved 
in targeting proteins to cell membranes. C2 domains are 
modular lipid-binding domains found in a variety of pro- 
teins with functions that include vesicular transport, 
GTPase regulation, lipid modification, protein phosphor- 
ylation, and ubiquitylation [75]. The most extensively 
studied proteins included in the last group are the 
synaptotagmins (SYT). There are reports describing the 
involvement of C2-domain proteins in plant responses 
to abiotic and biotic stresses. Yang et al. [26,76] showed 
that the Arabidopsis C2-domain protein, BAPl, acts as a 
negative regulator of programmed cell death induced by 
biotic stimuli. This gene is differentially expressed in the 
present microarray experiment and in the Yariv mic- 
roarray experiment described by Guan and Nothnagel 
[44]. The third most up-regulated gene in the present data 
set is a calcium-dependent lipid binding protein, CaLB 
(At3g57880), a C2 protein which is also up-regulated in the 
Yariv array. It seems relevant to highlight that the C2 
calcium-dependent membrane targeting protein Atlg70810 
was identified as an AGP6 interactor in the yeast- two hy- 
brid assay. Furthermore, in a yeast two-hybrid assay 



performed on Nicotiana alata [77] to identif)^ pollen pro- 
teins potentially involved in deciphering chemical signals 
provided by the pistil, two pistil-specific AGPs from the 
transmitting tract were used to hybridize against a pollen 
cDNA library. The authors found three pollen proteins 
that interacted with those pistil AGPs, one of which was a 
C2 domain-containing protein (NaPCCP). 

Analysis of the present results further identified a pos- 
sible involvement of AGP6 and AGPll in membrane traf- 
ficking, and that a plant retromer may be present in pollen 
tubes where the movement of vesicles is very important. A 
retromer is a multiprotein complex that is strongly con- 
served in eukaryotes. It is involved in the recycling of 
transmembrane receptors, which mediate the transport of 
vacuolar/lysosomal hydrolases [78]. In yeast, retromers 
consist of a large subunit with three vacuolar protein sort- 
ing (VPS) proteins, VPS35p, VPS29p and VPS26p. Homo- 
logs of these proteins are present in mammals and in 
plants. Recently, Oliviusson et al. [37] showed that the 
three VPS proteins, namely VPS35 (At3g51310), VPS29 
(At3g47810), and VPS26 (At5g53530) of the plant retro- 
mer complex are localized to multivesicular bodies (MVBs) 
in tobacco BY2 cells. MAGI, a homolog of VPS29, is dif- 
ferentially expressed in the present array, as is VPS26. 
Besides, VPS29 is also differentially expressed in the Yariv 
microarray data and MAGI is one of the AGPll interac- 
tors in the yeast-two hybrid experiment. Recent data 
strongly suggest that the retromer complex has conserved 
a function in mediating retrograde trafficking [38] and that 
VPS29 is involved in the cycling of certain PM proteins 
and required for the establishment of cell polarity during 
organ initiation in plants [66]. 

A working model for the role of AGPs in pollen tube 
development which integrates the data obtained in this 
work is proposed (Figure 6). 

Conclusions 

It is known that plant cells recycle pectins and AGPs 
[80]. Pectin cross-linking calcium may be removed and 
pectins recycled, maintaining the loosened walls essen- 
tial for growth. AGPs can also be recycled, to maintain 
their concentration level at the pollen tube apex neces- 
sary for them to perform a signaling role. 

All the data obtained in the present work emphasizes 
the remodeling of the plasma membrane via endocytosis. 
Endosomal trafficking pathways are emerging as funda- 
mental regulators of the wall physiology, involved in 
multiple signaling pathways and developmental pro- 
cesses. The interaction of AGP6 and AGPll with mem- 
bers of the pollen tube endosome machinery provides 
strong evidence for the recycling of these cell wall pro- 
teins (Figure 6). 

It is our understanding that as we continue to unravel 
the signal transduction processes involved in intercellular 
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Figure 6 A model for the role of AGPs in pollen tube growth. (A) Accepted model for pollen tube growth [70,79], illustrating vesicle 
exocytosis (1), smooth endocytosis (2), clathrin vesicle endocytosis (3), the pathway into degradation (4) and the recycling pathway (5). The apical 
dome (red), the exocytosis region (green), the multivesicular body (purple) and the vacuole (blue) are shown. (B) Proposed model for AGP 
signaling role in the pollen tube apical dome. AGPs act as receptors for extracellular signals and interact with transmembrane proteins, possibly 
receptor kinases or C2 domain-containing proteins. These interactions lead to opening of calcium channels, triggering various intracellular events. 
During pollen tube growth, AGPs are recycled by endocytosis, and either reused or sent for degradation, through multivesicular bodies. (C) 
Proposed model for the sub-apical zone. The hardening of cell wall pectins resulting from complexing with calcium hinders AGP contact with 
putative signal molecules, rendering AGPs inactive. 



communication during pollen/pistil interaction, we will 
discover more about the function of AGPs. 

Methods 

Plant material and growth conditions 

A Ds-tagged agp6 agpll homozygous double null mutant 
line [11] obtained from crossing RIKEN lines Ds54-4754-l 
and Dsll-4025-1 (RIKEN GSC Arabidopsis Ds transposon 
tag line collection, [81]), was used for pollen harvest and 
plants from the ecotype Nossen (No-0 Line N3081 NASC) 
were used as control. All plants were germinated and 
grown in half strength Murashige and Skoog (MS) 
medium complemented with 0.7% agar. Plantlets were 
then transferred to soil and kept in a growth chamber at 
22°C under long days (16 h light/8 h dark), irradiance of 
130 (imol m'^ s'^ and 60% relative humidity. 

Pollen isolation and pollen tube culture 

Pollen was collected from newly opened flowers and ger- 
minated according to Dardelle et al. [82]. Forty open 
flowers of Arabidopsis were introduced in each 1.5 mL 
microtube, and gently vortexed with 1 mL of pollen 
growth medium (5 mM KCl, 5 mM CaC^, 1 mM MgS04, 
0.01% H3BO3, 10% Sucrose, pH 7.4) for 5 min. The 
flowers were then removed with forceps, the microtubes 



centrifuged for 3 min at 3000 rpm and the supernatant 
discarded. The pelleted pollen destined for pollen tube 
culture was combined in groups of three and resuspended 
in 1.25 mL of fresh pollen growth medium. The prepar- 
ation was then transferred to 20 mm diameter flat bottom 
flasks and incubated at 22°C in the dark for 8 h. The ger- 
minated pollen tubes were then washed through a 60 (im 
mesh nylon net filters (Millipore NY6004700) with growth 
medium in order to discard ungerminated pollen grains. 
Pollen tubes from approximately 1600 flowers were used 
to obtain each independent replicate for agp6 agpll and 
for wild type. 

RNA isolation, target synthesis and hybridization to 
affymetrix GeneChips 

Total RNA was extracted using the RNeasy Mini Kit 
(Qiagen, Hilden, Germany). Concentration and purity was 
determined by spectrophotometry and integrity was con- 
firmed using an Agilent 2100 Bioanalyzer with a RNA 
6000 Nano Assay (Agilent Technologies, Palo Alto, CA). 

RNA was processed for use on Affymetrix (Santa Clara, 
CA, USA) GeneChip Arabidopsis ATHl Genome Arrays, 
according to the manufacturers GeneChip 3/ IVT Express 
kit user manual. Briefly, 100 ng of total RNA containing 
spiked in Poly-A RNA controls was used in a reverse 
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transcription reaction (GeneChip 3/ IVT Express Kit; 
Afiymetrix) to generate first-strand cDNA. After second- 
strand synthesis, double-stranded cDNA was used in a 
16 h in vitro transcription (IVT) reaction to generate 
aRNA (GeneChip 3/ IVT Express Kit; Aftymetrix). Size 
distribution of the aRNA and fragmented aRNA, respect- 
ively, was assessed using an Agilent 2100 Bioanalyzer with 
a RNA 6000 Nano Assay. 

15 (ig of fragmented aRNA was used in a 300 \A 
hybridization cocktail containing added hybridization con- 
trols. 200 (il of mixture was hybridized on arrays for 16 h 
at 45°C. Standard post hybridization wash and double-stain 
protocols (FS450_0004; GeneChip HWS kit, Affymetrix) 
were used on an Affymetrix GeneChip Fluidics Station 
450. Arrays were scanned on an Aftymetrix GeneChip 
scanner 3000 7 G. 

GeneChip data analysis 

Scanned arrays were analyzed first with Affymetrix 
Expression Console software to obtain Absent/ Present 
calls and to assure that all quality parameters were in the 
recommended range. Subsequent analysis was carried out 
with DNA-Chip Analyzer (dChip) 2010 (http://www. 
dchip.org, Cheng Li Lab, Harvard). The 6 arrays were nor- 
malized to a baseline array with median CEL intensity by 
applying an Invariant Set Normalization Method [83]. 



Normalized CEL intensities of the arrays were used to ob- 
tain model-based gene expression indices based on a PM 
(Perfect Match)-only model [84]. Replicate data (tripli- 
cates) for the conditions were weighted gene-wise by using 
inverse squared standard error as weights. Genes com- 
pared were considered to be differentially expressed if the 
90% lower confidence bound of the fold change between 
experiment and baseline was above 1.3. The lower confi- 
dence bound criterion means that we can be 90% 
confident that the fold change is a value between the lower 
confidence bound and a variable upper confidence bound. 
Li and Wong [83,84] have shown that the lower confi- 
dence bound (LBFC) is a conservative estimate of the fold 
change and therefore more reliable as a ranking statistic 
for changes in gene expression. 

Quantitative real-time PGR (qPCR) 

Samples from the same RNA preparations used in the 
microarray experiment were reverse transcribed using 
Promega Reverse Transcription System and poly(dT)i2_i8 
to prime the reactions. cDNA was amplified using the iQ™ 
SYBR® Green Supermix on the iQ™5 Real-Time PCR 
Detection System (Biorad). 

All qPCR reactions were run in duplicates. Thermocycle 
settings were as follows: Initial denaturation of 3 min at 
95°C, followed by forty cycles, each consisting of 10 s at 



Table 5 Oligonucleotide primer sequences used in qPCR and in reverse-transcription PCR assays 



AGI ID 


primer ID 


Forward primer 


Reverse primer 


At4g05320 


UBQIO 


gctccgacaccattgacaac 


acgcaggaccaagtgaagag 


Atlg71380 


ATCEL3 


tctccaagtcattgctcttcttcc 


cgttgtctcctgcgtcatagtac 


At3g60570 


ATEXPB5 


gcaaacggtgatgggaacttcg 


ggacacggcggaggtaagc 


AtSgllllO 


ATSPS2F 


tggtggtgttcgtgggagattcag 


ttagcctcggtgatgttgggactg 


Atlg53540 


HSP17.6C-CI 


ggcaaacgcacccgctatg 


ttcacctcttccttcctcagtcc 


AT2g46500 


AT2g46500 


agacggctcaaacgctcagaatc 


cggatagggcttcgaggaatgc 


At2g39890 


PROTl 


atggcgagaggcgggtac 


gtggttggctagaatgaatgtgag 


Atlg50490 


UBC20 


agatcctccggcgtctaatgg 


tgcttcctgttattgtccctttcc 


At5g60250 


AT5g60250 


cggttacagatggaggaggcactc 


aggacaccacacgtagccacaatc 


Atlg68610 


Atlg68610 


ctctaacgaccaaccaagccaag 


acatcgctccgctcacacc 


At4g27960 


UBC9 


gtttcaccaccctttcttc 


aaatcccacgatcaaattcc 


At2g 19970 


At2g 19970 


gagacattctgatggaccttacg 


tccgagacttaacgattgattgg 


Atlg65760 


Atlg65760 


cgatgattcctacattagcagac 


aagacgcaacgggtaacg 


At2g22340 


At2g22340 


gcgacggtgggatttcaggaactg 


gatgggagcggcaacggatgtg 


At5g28380 


At5g28380 


aagcacactgcgaccaaggc 


cagcggctcatggatctcactac 


At4g09950 


At4g09950 


ggaggatgtgaaggagcaattagc 


tcttgttgagttcggtgcgtaac 


Atlg 12840 


DET3 


cggcgttcttggcatgtgtc 


agcaaggttgatagtgaaggagac 


Atlg51490 


BGLU36 


gccgtactctctcgctgtcaaag 


atgagccagaagttcgtgatgtcc 


At3g57880 


C2 domain-containing protein 


caggatgaggtatrgacaggttgag 


gcggcaatcaagcagaacaag 


At4g20780 


CML42 


cccaagcctaaacgcacttcg 


acggtggatttgagatcggagag 


At2g 19970 


CAP (Antigen 5) 


gagacattctgatggaccttacg 


tccgagacttaacgattgattgg 


At2g 19980 


CAP (allergen V5) 


gcacagaggtacgctaacg 


ggtggcataattgtaataaggc 
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Table 6 Oligonucleotide primer sequences used in the DNA constructions for the Y2H assays 



AGI ID 


primer ID 


Forward primer 


Reverse primer 


Atbg 14380 


AGP6 


ggggacaagtttgtacaaaaaagcaggctcggccgacgctccctcagcttctcc 


ggggaccactttgtacaagaaagctgggtcctaactcttgggtgactctgcagtgg 




AGPl 1 


y y y y dLddy LLLLj IdLddddddy Ldy y LLLy y LLy d ly LdLLLLLdy Liy LdLL 


y y y y dLLdLLLLy idLddydddyLiy yy iLLLddLLLLiyyy lydLiLyyLyyL 


Atlg26630 


FBR12 


ggggacaagtttgtacaaaaaagcaggctcgaccttcctcttcccctcc 


ggggaccactttgtacaagaaagctgggtcctagtgcaaatagcaatgaatatc 


At4g05320 


UBQIO 


ggggacaagtttgtacaaaaaagcaggctcgaaatcttaaaaactttctctc 


ggggaccactttgtacaagaaagctgggtcctaaaacaaaagaagcacagataat 


AtSg 18040 


MAPK9 


ggggacaagtttgtacaaaaaagcaggctcggcgaaaagtttctcccttg 


ggggaccactttgtacaagaaagctgggtcctatgaagaaaacacacactttaac 


Atlg70810 


CaLB 
domain 
family 
protein 


ggggacaagtttgtacaaaaaagcaggctcgaaaagagtcagagccgc 


ggggaccactttgtacaagaaagctgggtcctaaaaccaaaacgatttagg 








At2g30020 


AP2C1 


ggggacaagtttgtacaaaaaagcaggctcgcgaaatcgaaatcaaaatatag 


ggggaccactttgtacaagaaagctgggtcctaaatttcggccaatgctcg 



95°C, 30 s at 56°C and 30 s at 72°C. After each run, a dis- 
sociation curve was acquired to check for amplification 
specificity by heating the samples from 60 to 95°C. 

Serial dilutions of both agp6 agpll and wild type 
pollen tube cDNA were used to determine the efficiency 
curve of each primer pair. The primers used are listed in 
Table 5. UBQIO was used as internal reference gene. At 
the end of the PGR cycles, the data was analyzed with 
the iQ5 2.0, Standard Edition Optical System Software 
v2.0. 148.060623 (Biorad), using the Livak calculation 
method [85]. 

Semi quantitative reverse transcription (RT)-PCR 

cDNA samples used in qPCR experiments were also 
used in the RT-PCR assays. RT-PCR reactions were done 
on 8 genes, and two reference genes, UBC9 and UBQIO 
(Table 5). PGR reactions were set up with DreamTaq 
DNA Polymerase (Fermentas). Each PGR reaction was 
initiated with a 5 min denaturation at 95°G followed 
cycles of 30 s at 95°G, 30 s at 56°G (for all primer pairs) 
and 30 s at 72°G. Samples were taken after 25, 30 and 
35 cycles. cDNA quantities were normalized for each 
sample using the two reference genes. Band intensity 
was assessed using Kodak DG120 Gel Electrophoresis 
Analysis System. 

Yeast Two-hybrid library screening assays 

The predicted peptide cores of AGP6 and AGPll, fol- 
lowed by a stop codon, were individually cloned into the 
pGBTKT? vector (Glontech, Palo Alto, GA), modified at 
the Xmal site to include a Gateway cassette. The bait 
constructs were introduced in the a-Y187 strain [86]. 
YPAD, SD, and appropriate dropout media have been 
described previously [87]. 

The bait strains were mated with a normalized Arabidopsis 
thaliana total plant cDNA library cloned in pGADrec 
(Sommer and Masiero, unpublished data) and introduced 
in the a-AH109 yeast strain (Glontech). 

Golonies that grew on all selective media (-Trp-Leu- 
Adenine-His and supplemented with 15 mM 3-Amino- 



1,2,4-triazole) were further characterized. The pGADrec 
plasmids were rescued and cDNA insert was amplified by 
PGR and sequenced. 

Yeast Two-hybrid interaction assays 

cDNA of each candidate gene was used to perform con- 
firmation of the interaction. RNA was extracted from 
Arabidopsis inflorescences containing flowers in several 
stages of development using the RNeasy Mini Kit (Qiagen). 
cDNA was synthesized with the RevertAid First Strand 
cDNA Synthesis Kit (Thermo Scientific, Waltham, MA, 
USA). Full-length cDNA of these genes was amplified by 
PGR (Table 6) with Phusion High-Fidelity DNA Polymer- 
ase (New England Biolabs, Ipswich, MA, USA) and puri- 
fied from agarose gel using the GeneJET Gel Extraction 
Kit (Thermo Scientific). In case of alternative splicing 
forms, the splice selected was the one corresponding to 
the mature form according to TAIR. The cDNA sequence 
of each gene was individually cloned into the pGADT7 
vector (Glontech). 

Each bait/prey pair was introduced in the a- AH 109 
yeast strain (Glontech), and as a control for autoactiva- 
tion false-positives, each bait was also co-transformed 
into the yeast strain with the empty AD vector, and each 
prey was co-transformed with the empty BD vector. 
Bait/prey pair colonies that grew on all selective media 
(-Trp-Leu-Adenine-His and supplemented with increas- 
ing concentrations of 1 mM to 2.5 mM 3-Amino- 1,2,4- 
triazole) were considered positive for interaction. 

Additional files 



Additional file 1: agp6 agpll mutant microarray data. 

Additional file 2: Differentially expressed genes common to agp6 
agpl 1 mutant and treatment with Yariv phenylglycoside reagent 
for 1 h and 10 h. 
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